{
"cells": [
{
"cell_type": "markdown",
"id": "dd3d8d19",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"# LOv rewriting rules in Perceval\n",
"\n",
"The aim of this notebook is to rewrite a circuit using rewriting rules based on the article: *LOv-Calculus: A Graphical Language for Linear Optical Quantum Circuits*.\n",
"\n",
"We show how to use these rewriting rules to generate unique triangular normal forms."
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "96bdb88a",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"import perceval as pcvl\n",
"import perceval.lib.phys as phys\n",
"from perceval.algorithm.optimize import optimize\n",
"from perceval.algorithm.norm import frobenius\n",
"import random"
]
},
{
"cell_type": "markdown",
"id": "ba837e96",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"This is the first rewrite rule used in this noteobok. It is the rule 37 in the article."
]
},
{
"cell_type": "markdown",
"id": "09923a89",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
""
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "ede320ce",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pattern1=pcvl.Circuit(3, name=\"pattern1\")//(0,phys.PS(pcvl.P(\"phi0\")))//(0,phys.BS(theta=pcvl.P(\"theta1\")))//(0,phys.PS(pcvl.P(\"phi2\")))//(1,phys.PS(pcvl.P(\"phi1\")))//(1,phys.BS(theta=pcvl.P(\"theta2\")))//(0,phys.BS(theta=pcvl.P(\"theta3\")))\n",
"pattern1._color = \"lightgreen\"\n",
"pcvl.pdisplay(pcvl.Circuit(3).add(0,pattern1,False), recursive=True)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "1e1ee73a",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rewrite1=pcvl.Circuit(3, name=\"rewrite\")//(0,phys.PS(pcvl.P(\"beta2\")))//(1,phys.PS(pcvl.P(\"beta1\")))//(1,phys.BS(theta=pcvl.P(\"alpha1\")))//(0,phys.BS(theta=pcvl.P(\"alpha2\")))//(1,phys.PS(pcvl.P(\"beta3\")))//(1,phys.BS(theta=pcvl.P(\"alpha3\")))//(0,phys.PS(pcvl.P(\"beta4\")))//(1,phys.PS(pcvl.P(\"beta5\")))//(2,phys.PS(pcvl.P(\"beta6\")))\n",
"rewrite1._color = \"lightgreen\"\n",
"pcvl.pdisplay(pcvl.Circuit(3).add(0,rewrite1,False), recursive=True)"
]
},
{
"cell_type": "markdown",
"id": "d4197839",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"Let us implement now the rule number 1."
]
},
{
"cell_type": "markdown",
"id": "171cff3a",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
""
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "5468fdb9",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"pattern2=pcvl.Circuit(1, name=\"pattern\")//phys.PS(pcvl.P(\"phi1\"))//phys.PS(pcvl.P(\"phi2\"))\n",
"rewrite2=pcvl.Circuit(1, name=\"rewrite\")//phys.PS(pcvl.P(\"phi\"))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "e428775f",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pcvl.pdisplay(pcvl.Circuit(1).add(0,pattern2,False), recursive=True)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "e2bdd9f4",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pcvl.pdisplay(pcvl.Circuit(1).add(0,rewrite2,False), recursive=True)"
]
},
{
"cell_type": "markdown",
"id": "6ceae57c",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"The third rule used in this notebook is the following one:"
]
},
{
"cell_type": "markdown",
"id": "246fc6ce",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
""
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "d056ab6c",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pattern3=pcvl.Circuit(2, name=\"pattern3\")//(1,phys.PS(pcvl.P(\"phip\")))//(0,phys.BS(theta=pcvl.P(\"theta\")))\n",
"pattern3._color = \"pink\"\n",
"pcvl.pdisplay(pcvl.Circuit(2).add(0,pattern3,False), recursive=True)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "6ef8057d",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rewrite3=pcvl.Circuit(2, name=\"rewrite3\")//(0,phys.PS(pcvl.P(\"phi1\")))//(0,phys.BS(theta=pcvl.P(\"theta\")))//(0,phys.PS(pcvl.P(\"phi2\")))//(1,phys.PS(pcvl.P(\"phi3\")))\n",
"rewrite3._color = \"pink\"\n",
"pcvl.pdisplay(pcvl.Circuit(2).add(0,rewrite3,False), recursive=True)"
]
},
{
"cell_type": "markdown",
"id": "6c7efe31",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"And the fourth rule is the rule 38 in the article."
]
},
{
"cell_type": "markdown",
"id": "4cae9637",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
""
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "597bc933",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pattern4=pcvl.Circuit(2, name=\"pattern4\")//(0,phys.BS(theta=pcvl.P(\"theta1\")))//(0,phys.PS(pcvl.P(\"phi1\")))//(0,phys.BS(theta=pcvl.P(\"theta2\")))\n",
"pattern4._color = \"orange\"\n",
"pcvl.pdisplay(pcvl.Circuit(2).add(0,pattern4,False), recursive=True)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "1866e2c5",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rewrite4=pcvl.Circuit(2, name=\"rewrite4\")//(0,phys.PS(pcvl.P(\"beta1\")))//(0,phys.BS(theta=pcvl.P(\"alpha1\")))//(0,phys.PS(pcvl.P(\"beta2\")))//(1,phys.PS(pcvl.P(\"beta3\")))\n",
"rewrite4._color = \"orange\"\n",
"pcvl.pdisplay(pcvl.Circuit(2).add(0,rewrite4,False), recursive=True)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "deafebff",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"a=pcvl.Circuit.generic_interferometer(4, lambda idx:pcvl.Circuit(2)//phys.PS(phi=random.random())//phys.BS(theta=random.random()), depth=8, shape=\"rectangle\")\n",
"pcvl.pdisplay(a, recursive=True, render_size=0.7)"
]
},
{
"cell_type": "markdown",
"id": "f13dd7d7",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"## Normalizing Circuit"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "0f5ee748",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"import drawSvg as draw\n",
"import copy"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "df109c0f",
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"matching pattern pattern\n",
"matching pattern pattern3\n",
"matching pattern pattern4\n",
"matching pattern pattern1\n",
"matching pattern pattern\n",
"matching pattern pattern3\n",
"matching pattern pattern4\n"
]
}
],
"source": [
"reverse = []\n",
"direct=[]\n",
"def draw_frame(a):\n",
" if isinstance(a, pcvl.Circuit):\n",
" d = pcvl.pdisplay(a, recursive=True, render_size=0.6)\n",
" reverse.insert(0, d)\n",
" direct.append(d)\n",
" return d\n",
" return a\n",
"\n",
"rules = [(pattern1, rewrite1, \"lightgreen\"), (pattern2, rewrite2, \"lightblue\"),\n",
" (pattern3, rewrite3, \"pink\"), (pattern4, rewrite4, \"orange\")]\n",
"\n",
"with draw.animate_jupyter(draw_frame, delay=0.1) as anim:\n",
" anim.draw_frame(a)\n",
" while True:\n",
" found = False\n",
" for pattern, rewrite, color in rules:\n",
" start_pos = 0\n",
" while True:\n",
" print(\"matching pattern\", pattern._name)\n",
" matched = a.match(pattern, browse=True, pos=start_pos)\n",
" if matched is None:\n",
" break\n",
" print(\"matching ok\", matched.v_map)\n",
" idx = a.isolate(list(matched.pos_map.keys()), color=color)\n",
" anim.draw_frame(a)\n",
" for k, v in matched.v_map.items():\n",
" pattern[k].set_value(v)\n",
" v = pattern.compute_unitary(False)\n",
" print(\"optimizing rewrite\",rewrite._name)\n",
" res = optimize(rewrite, v, frobenius, sign=-1)\n",
" print(\"found params with distance\", res.fun)\n",
" subc = rewrite.copy()\n",
" found = True\n",
" a.replace(idx, subc, merge=False)\n",
" anim.draw_frame(a)\n",
" a.replace(idx, subc, merge=True)\n",
" pattern.reset_parameters()\n",
" rewrite.reset_parameters()\n",
" anim.draw_frame(a)\n",
" start_pos = idx\n",
" if not found:\n",
" break"
]
},
{
"cell_type": "markdown",
"id": "954dfac2",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"This representation is exactly the normal form that we wanted to obtain !"
]
},
{
"cell_type": "markdown",
"source": [
"## Reference\n",
"\n",
"A. Clément, N. Heurtel, S. Mansfield, S. Perdrix, B. Valiron. LOv-Calculus: A Graphical Language for Linear Optical Quantum Circuits, [arXiv:2204.11787](https://arxiv.org/abs/2204.11787) (2022)."
],
"metadata": {
"collapsed": false,
"pycharm": {
"name": "#%% md\n"
}
}
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.4"
}
},
"nbformat": 4,
"nbformat_minor": 5
}